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^ ■ Abstract 

m 

Measurements of heavy quark production at a possible future ep facility, THERA, would 
. provide valuable information on the structure of the photon. QCD cross-section predictions are 

[ made at LO for both charm and beauty production and their sensitivity to current parametri- 

' sations of the photon parton densities is investigated. 

m ■ 

1 Introduction 

r~| ' Heavy quark production can provide information on the structure of the real photon as well 

as provide stringent tests of QCD by itself. In this paper the potential of the possible future 
Q^l ep collider, THERA, to give new information on the real photon structure through charm and 

beauty quark photoproduction is studied. The potential of the collider is also compared with 
other existing or planned colliders. 

Photoproduction processes are dominated by the exchange of quasi-real photons (virtuality 
X ■ of the photon ^ iQeV^). Their energy distribution can be described by an effective (real) 

. photon energy spectrum, i.e. using the Weizsacker- Williams approximation. A hard scale is 

provided by a large transverse momentum or the mass of the produced heavy quark Q. In this 
paper the region pT > mq is studied. In LO QCD, two processes contribute to the photopro- 
duction of heavy quarks: the direct process, in which the photon couples directly to a parton 
in the proton, and the resolved process, in which the photon acts as a source of partons, one of 
which scatters from a parton in the proton. Heavy quarks can be produced in the final state 
or they already exist in the initial state as partons in the photon and proton, depending on the 
scheme used in the calculations. 

If not stated otherwise the computations were performed with the THERA electron and 
proton beams energies of 250 GeV and 920 GeV, respectively. In the first section results are 
presented on LO QCD calculations of inclusive charm quark production, ep — >■ c/cX, at the 
THERA collider. In particular the possible information this process can reveal on the structure 
of the photon is discussed. The second section uses MC expectations to investigate events where 
at least one heavy quark and two jets are produced. Further requirements on a final state muon 
from a semi-leptonic decay are also imposed and sensitivity to current parametrisations of the 
photon investigated. 
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2 Leading order QCD predictions of charm quark production 



LO calculations were performed in two schemes which differ by the number of quark flavours 
which are considered to be part of the structure of the photon and proton, and can thus take 
part in the process as partons. The massive scheme, the so-called Fixed Flavour Number Scheme 
(FFNS), assumes that both the photon and proton consist of gluons and only light {u, d and 
s) quarks whilst heavy quarks are produced in the final state. The massless Variable Flavour 
Number Scheme (VFNS), considers heavy quarks as active flavours. Therefore such partonic 
processes, Q + g^Q + govQ + Q^Q + Q, with Q denoting a heavy quark, are also allowed. 
(Events in which the heavy quark is one of the active partons are also referred to as "flavour- 
excitation".) In the FFNS scheme the mass of the produced charm quark is always kept nonzero 
in the calculations while in the VFNS all flavours are treated as massless. 

In the massive (FFNS) calculations the number of active flavours {Nf) is taken to be 3. In the 
massless (VFNS) scheme it varies from 3 to 4 depending on the value of the hard (factorisation 
and renormalisation) scale, /i. The charm quark is included in the computation provided that 
^ > TTT-c with the mass of the charm quark set to rUc = l.SGeV. Also the QCD energy scale, 
Aqcd, which appears in the one loop formula for the strong coupling constant is affected by 
the change of the number of active flavours, so is denoted as Aq^^. The scale is taken to be 
A^fjjj = 232 MeV and A'^^D = 200 MeV as in the GRVQ parametrisation of the photon ||]. 
If not stated otherwise is taken to be the transverse mass of the produced charm quark. 



rriT = yiTi'l + Pr-: ^-^d the differential cross-section, (P a / drjclp"^ , with pr = 10 GeV is calculated. 
This means that ji ^ pT and therefore Nj = 4 in case of the VFNS calculations. The spectrum 
of the quasi-real photons taking part in the process is described by the Weizsacker- Williams 
(EPA) formula. 



A/e(y) = £ 



where Q^ax — 4GeV , Q^in = "T-e?/ /(I ~ v) ^"^^ U is the fraction of the electron's four- 
momentum carried by the photon {y = E^/Eg). No requirements were placed on y. Finally 
the LO CTEQ5L Q and GRV were used as the parton parametrisations of the proton and 
photon, respectively. 

To show the power of the THERA collider for heavy quark physics it is useful to compare 
the cross-section for charm production that can be achieved in THERA with cross-sections 
expected for other existing or planned accelerators. In Fig. [l|, a comparison with results from 
the ep HERA {Ee = 27.5 GeV, Ep = 920 GeV), LEP {Ee = 90 GeV), a Linear Colhder (LC) and 
a Photon Collider (PC) both with Eg = 250 GeV are presented for two options of THERA, ep 
and 7p. Real photons in the PC and jp THERA option can be achieved through the Compton 
backscattering laser photons from the electron beam 1^, ^ (the full theoretical spectrum for 
< y < 0.83 was used). It can be seen that THERA provides the largest cross-sections over 
the whole range in r] and has the potential to probe regions kinematically inaccessible at other 
colliders. 

It has been demonstrated elsewhere |Q, that there is a large dependence of the cross-section 
for heavy quark production on the choice of scheme for the case of HERA energies. It is also 
true for the THERA collider. Direct and resolved contributions to the cross-section for charm 



^The old GRV{'92) parametrisation was used in this analysis. It treats heavy quarks as massless above the 
heavy quark threshold region {W 3> Smg), and moreover latest experimental results Q indicate that its gluon 
distribution is closest to the data. 
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Figure 1: Inclusive charm quark production in processes including real photon interactions 
calculated in the VFNS scheme for pT = 10 GeV for five accelerators: THERA (ep and jp 
options) and HERA {Ee = 27.5 GeV, Ep = 920 GeV) {ep c/cX) , LEP {Ee = 90 GeV) & LC 
{Ee = 250 GeV) (e+e" — > e~^e~c/cX) with quasi-real WW photons and PC (e+e~ e^e^c/cX) 
with real photons achieved through the Compton backscattering laser photons from the electron 
beam and Ef. = 250 GeV |Q, |5|). CTEQ5L proton and GRV LO photon parton parametrisations 
were used. 

production obtained using the two discussed schemes clearly show this dependency, as displayed 
in Fig. |2|. Note that the resolved and direct cross-sections are similar for the VFNS scheme 
whereas in the FFNS scheme, the resolved cross-section is small compared to that of the direct 
contribution. The small difference between the direct cross-sections for the two schemes is 
due to the mass term in the matrix elements and phase space integral. Figure ^ presents the 
cross-sections with their dependence on the choice of the ^ varied from l/2mr to 2mT- The 
dependence on the choice of scheme is most significant for resolved processes, varying by up to 
several orders of magnitude in the cross-section. In the case of direct processes the variation is 
relatively small and the subsequent total cross-section prediction variation is therefore dominated 
by the difference between the two schemes for the resolved process. The variation in the choice 
of scale leads to small effects in the cross-section in comparison with the scheme dependence as 
it only influences the value of and the parton densities and not the set of allowed partonic 
processes. 

To test the sensitivity of inclusive charm production to the form of the parton densities 
in the photon, cross-sections for this process were computed with three different LO parton 
parametrisations of the photon: GRV GS96 0] and SaSlD [^. Sensitivity to the form 
of the current parametrisations of the photon structure function at the THERA collider is 
larger than at HERA as shown in Fig. ^ (the results of the VFNS calculation). Note that 
the SaSlD and GRV parametrisations give similar results. Figure |^ shows the uncertainty 
resulting from the variation of the hard scale by a factor of two compared to the central value. 
Here it is demonstrated that the differences between the cross-sections obtained using three 
parton parametrisations of the photon are significant for GRV (or SaSlD) and GS96 in the 
region —2.5 < < 0, despite this large theoretical uncertainty. 
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Figure 2: THERA. (a) Direct and (b) resolved contributions to the differential cross-section 
of ep — > c/cX production calculated for pT = 10 GeV in the VFNS (solid lines) and FFNS 
(dashed lines) schemes using the CTEQ5L and GRV LO parton parametrisations of the proton 
and photon, respectively. 




Figure 3: THERA. Differential cross-section of ep — > c/c X production calculated for px = 10 
GeV in the VFNS (solid lines) and FFNS (dashed lines) schemes with three different values of 
the hard scale; n = 1/2, 1 and 2 x rriT, using the CTEQ5L and GRV LO parton parametrisations 
of the proton and photon, respectively. 
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The differences between two current parton parametrisations of the photon, GRV and SaSlD, 
have previously been compared for LC and PC colhders More specificahy, the relative 
difference '^^^grv^^^ of the cross-section, cPa/dridpj^{ep — > c/cX), was calculated in both 
VFNS and FFNS schemes. Values of ~ 0.1 — 0.4 were found with the PC having the larger 
value. Similar numbers (up to 0.3) to those for the LC and PC could be achieved at THERA 
in the forward 7] region. The massive FFNS gives the same results for the considered sensitivity 
as those of the massless VFNS scheme. 




Figure 4: Differential cross-section for ep — > c/cX for pT = 10 GeV calculated with the CTEQ5L 
parton parametrisation of the proton and three various LO parton parametrisations of the photon 
in the massless (VFNS) scheme. A) Lower lines correspond to the HERA collider energies 
(27.5x920 GeV), while upper to the THERA collider energies. B) Lower and upper lines for 
each photon parametrisation show respectively the lowest and greatest values of the cross-section 
calculated with different values of the hard scale /x = l/2,l, &2x itlt in the GRV case and 
|U = 1, 2 X rriT in the GS96 case. 

Sensitivity to the photon parton densities can be also tested through the ratio of the resolved 
to the direct photon contributions. This is shown in Fig. ^ for d?a/dr]dp'^ in comparison with 
the results for the HERA collider. Again the results for THERA show a larger resolved cross- 
section and increased sensitivity to the form of the current parametrisations. The sensitivity is 
also large in the central ij region where there is, experimentally, a large acceptance. 

3 MC predictions of dijet production with charm and beauty 
quarks 

The dijet system can be formed from the production of a heavy quark pair such as in the boson- 
gluon fusion process. When the heavy quark is an active parton in the photon then the dijet 
system is formed from one heavy quark and a gluon. The other heavy quark in the event is 
produced at low transverse energy and forms part of the photon remnant. 

A natural variable to describe dijet production is the observable x^^, which is the fraction 
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Figure 5: The ratio of the resolved to the direct photon contribution of e'^e" e^e~ c/c X 
(for d?a/dr]dpj^) for pT=10 GeV calculated in LO with the CTEQ5L parton parametrisation of 
the proton and three various LO parton parametrisations of the photon in the massless (VFNS) 
scheme. A) HERA, B) THERA cohider. 



of the photon's energy producing the two jets of highest transverse energy [10 



' ~ 2yE, ' 



where rf^^ is the pseudorapidity of the jet and is the electron energy and y the fraction of 
the electron's energy carried by the photon. 

Due to the increase in the beam energy, there is a corresponding order of magnitude exten- 
sion in the minimum x°^^. This is shown in Fig. ^ in which two THERA scenarios (with the 
electron beam energy = 250 GeV and E^ = 400 GeV) are compared with HERA. A greatly 
increased cross-section possible at THERA over that at HERA for low— x"^*^, particularly in 
beauty production, is seen. Realistic scenarios for the cuts on the jet quantities have been made 
of transverse energy E";^*^'^ > 10,8 GeV and pseudorapidity if^^ < 2.5. The distribution is 
peaked at high-x"''^, consistent with direct photon processes, but has a significant cross-section 
at low-x"*^*^ arising from resolved photon events. 

Charm and beauty production at a 7p collider has also been considered and compared with 
the nominal ep scenario for THERA. The photon beam energy is assumed to be 0.8Ee, i.e. 
200 GeV assuming an electron beam of energy 250 GeV. The predictions for the ep and jp 
options with the same jet cuts are shown in Fig. |^. The cross-section prediction for the jp 
option is a factor of ~ 25 higher than the nominal ep option. This factor arises simply from the 
Weizsacker- Williams flux associated with the electron. 

Further requirements can be made such that the kinematic range closely corresponds to 
what could be measured. Therefore the observation of a final state muon from the semi-leptonic 
decay of a beauty quark is imposed. The momentum and angular requirements imposed on 
the muon are typical of those at HERA. Considering jets of slightly lower transverse energy 
(£';j?*^'^ > 7,6 GeV), the effect on the low-x°^*^ component of pseudorapidity cut is shown in 
Fig. p. It can be seen that for a pseudorapidity cut of rf"^^ > —1, the direct peak at high-x"*^*^ 
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Figure 6: The differential cross-section, da /dlogiox'^^ for (a) charm and (b) beauty in dijet 
photoproduction at HERA and THERA as predicted by the Herwig MC. The dashed and 
soHd hues show THERA with an electron energy of 400 GeV and 250 GeV respectively and the 
dot-dashed lines indicate the expectation from HERA. The kinematic range is the same in all 
three cases, namely, Q'^ < I GeV^, with two jets, > 10,8 GeV and rf''^ < 2.5. 




Figure 7: The differential cross-section, da /dlogiox'^'^ for (a) charm and (b) beauty in dijet 
photoproduction at for 7p and ep options at THERA as predicted by the Herwig MC. The 
solid lines show THERA with an electron energy of 250 GeV and the dashed lines the expectation 
with an photon energy of 200 GeV. The kinematic range is the same in all three cases, namely, 
two jets, Ei^^^'"^ > 10,8 GeV and rji"^ < 2.5. 
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is drastically reduced with respect to a cut of rf^^ > —2. This demonstrates that by applying a 
cut of r/J'^* > — 1 interactions of almost exclusively the resolved photon are being probed. 

This rather tight restriction on the pseudorapidity of the jets of — 1 < rf^^ < 2 is imposed 
and the sensitivity of the cross-section to the current parametrisations of the photon structure 
function investigated for THERA. In Fig. |^ predictions are shown for events where a beauty 
quark decays into a muon. Predictions are shown using both Herwig |11] and Pythia [^] MC 
generators for four different LO photon PDFs. At values of x°^'^=0.1, differences between the 
four parametrisations of up to 40% are seen. The predictions from GRV and GS96 generally 
have the largest cross-sections at low-x"'^^ and also rise more quickly. The large difference in 
the Pythia and Herwig MC predictions comes mainly from the different default treatments 
of Os , different scales and hadronisation and the use of massive matrix elements in Herwig and 
massless in Pythia for the generation of flavour-excitation processes. 




Figure 8: The differential cross-sections, (a) da/dx'^^ and (b) du/tilogiox"^'^ for beauty in dijet 
photoproduction at THERA as predicted by the Herwig MC. The events are selected in the 
range, <1 GeV^, 0.2 < y < 0.85 with two jets of transverse energy, E';^*^'^ > 7, 6 GeV. The 
points represent when the jets are restricted to — 1 < rf^^ < 2 and the solid lines when the jets 
are restricted to —2 < rf^^ < 2. A muon in the final state has also been required with transverse 
momentum, pj,>2 GeV and |r/^'[ < 2. 



4 Summary 

It has been shown that the photoproduction of heavy quarks at THERA has a much increased 
cross-section relative to HERA. With the expected luminosity of about 100 pb~^/year stringent 
tests of pQCD could be performed. It has also been shown that the cross-sections are sensi- 
tive to the structure of the photon, with the current parametrisations differing by significant 
amounts. The THERA collider also compares favourably with other current and planned ma- 
chines studying heavy quark production, with the 'jp mode for THERA providing the largest 
cross-section. 
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Figure 9: The differential cross-sections, (a) da/dx'^^ and (b) da / dlogiQx'^'^ for beauty in dijet 
photoproduction at THERA as predicted by the Herwig MC. The differential cross-sections, 
(c) da /dx'^^ and (d) da / dlogi^x'^^ for beauty in dijet photoproduction at THERA as predicted 
by the Pythia MC. Predictions for four different PDFs are shown; GS96-LO (open circles), 
GRV-LO (solid circles), SaS-lD (stars) and SaS-2D (diamonds). 



9 



Acknowledgements 

Some of use (PJ and MK) have been partly supported by the Pohsh State Committee for Scien- 
tific Research (grant 2P03B05119, 2000-2001). MK is also grateful to the European Commission 
50th Framework Contract HPRN-CT-2000-00149 and DESY for financial support. 



References 

[1] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 46 (1992) 1973. 
[2] C. Adloffat al., Euro. Phys. J. C 10 (1999) 363. 
[3] H. L. Lai et al., Euro. Phys. J. C 12 (2000) 375. 

[4] I. F. Ginzburg, G. L. Kotkin, V. G. Serbo and V. I. Telnov, Nucl. Instr. and Meth. 205 
(1983) 47; 

I. F. Ginzburg, G. L. Kotkin, S. L. Panfil, V. G. Serbo and V. I. Telnov, Nucl. Instr. and 
Meth. 219 (1984) 5. 

[5] V. I. Telnov, Nucl. Instr. and Meth. A 294 (1990) 72. 

[6] B. A. Kniehl, M. Kramer, G. Kramer and M. Spira, Phys. Lett. B 356 (1995) 539. 
[7] L. A. Gordon and J. K. Storrow, Nucl. Phys. B489 (1997) 405. 
[8] G. A. Schuler and T. Sj5strand, Z. Phys. C 68 (1995) 607. 

[9] P. Jankowki, M. Krawczyk and A. De Roeck, To be published in the proceedings of Work- 
shop on the Development of Future Linear Electron-Positron Colliders for Particle Physics 
Studies and for Research Using Free Electron Lasers, Lund, Sweden, 23-26 Sep 1999, Rolf 
Heuer, Francois Richard and Peter Zerwas Eds., hep-ph/0002169| ; and in preparation. 

[10] ZEUS Collaboration, M. Derrick et al., Phys. Lett. B 348 665 (1995) 

[11] G. Marchesini et al., Comp. Phys. Commun. 67 (1992) 465. 

[12] H.-U. Bengtsson and T. Sjostrand, Comp. Phys. Commun. 46 (1987) 43. 



10 



